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SUMMARY 


Computational  studies  on  the  jet  from  a forty-two  degree 
conical  liner  are  reported.  Three  different  modes  of 
initiation  were  studied  with  a view  to  predicting  the  best 
jet.  The  explosive  parameters  were  varied  systematically 
with  a view  to  determining  the  effects  of  each  parameter 
on  the  jet  formed. 
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EXPLOSIVE  OPTIMIZATION  FOR  A FORTY-TWO  DEGREE 
CONICAL  SHAPED  CHARGE 


INTRODUCTION 


In  optimizing  the  performance  of  a forty-two  degree  conical  shaped 
charge,  two  big  factors  are  the  explosive  chosen  and  the  mode  of 
initiation  for  the  explosive.  The  effects  of  both  of  these  variables 
can  be  studied  computationally  and  that  study  forms  the  subject  matter 
of  this  report . 

The  baseline  design  considered  for  this  study  was  the  forty-two 
degree  cone  shown  in  Figure  1.  (All  dimensions  are  in  centimeters). 

The  liner  is  made  of  copper  and,  for  computational  purposes,  elastic- 
plastic  effects  are  ignored  and  the  purely  hydrodynamic  Tillotson 
equation-of-state  is  used1 . The  explosive  is  pentolite  and  the 
reaction  products  are  treated  by  a gamma  law  equation-of-state.  The 
explosive  column  is  short,  being  only  about  3 centimeters  (O.M  cone 
diameters)  long. 

The  computations  were  carried  out  using  the  (two  dimensional) 

BRLSC  code  developed  by  Science,  Systems  and  Software  Corporation 
from  their  earlier  HELP  code  for  the  specific  purpose  of  computing 
shaped  charges1. 

In  Part  I,  three  different  modes  of  initiation  are  studied: 

A point  detonation  at  the  point,  P,  on  the  symmetry  axis  (see  Figure  1); 
a plane  wave  detonation  in  the  plane  of  the  points,  P and  Q,  (end 
face  of  the  explosive  column),  and  a circumferential  Initiation  in 
the  ring  obtained  by  rotating  Q around  the  symmetry  axis. 

In  Part  II,  three  different  hypothetical  explosives  were  used. 

The  baseline  calculation  was  the  plane  wave  initiated  cone  using 
pentolite  as  the  explosive.  By  varying  the  value  of  gamma  by  10^ 
a new  (hypothetical)  explosive  was  obtained.  The  results  of  this 
calculation  are  compared  to  the  baseline  calculation.  Similarly, 
two  more  explosives  are  obtained  by  varying  the  detonation  energy 
and  the  loading  density  by  10%.  While  these  explosives  correspond 
to  no  real  explosive,  they  each  serve  to  indicate  the  effect  of  the 
varied  parameter  on  the  jet  produced. 

PART  I - EFFECTS  OF  DIFFERENT  MODES  OF  INITIATION 


As  mentioned  above,  three  different  modes  of  initiation  are 
considered — point  Initiation,  plane  wave  initiation,  and 
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circumferential  initiation.  Tn  each  case  a detonation  wave  is 
formed  and  travels  down  the  explosive  until  it  collides  with  the 
liner.  Upon  collision,  a reflected  shock  arises  on  the  locus  of 
the  intersection  of  the  detonation  wave  and  the  liner.  At  10  psec 
after  initiation/the  pressures  in  the  explosive  are  as  shown  in  the 
pressure  maps  of  Figures  2,  3,  and  h . These  maps  are  based  on  an 
Eulerian  grid  which  divides  the  space  of  interest  into  fixed  cells. 

For  each  of  these  cells  a pressure  is  calculated.  In  order  to  present 
the  results  visually, a letter  is  assigned  to  various  pressure 
ranges  (see  Figures  2,  3,  and  4)  and  the  letter  appearing  in  the 
cell  indicates  the  range  to  which  the  cell  pressure  belongs.  The 
pressure  ranges  corresponding  to  each  letter  are  shown  at  the  bottom 
of  the  figure.  (The  ranges  are  not  the  same  for  all  three  figures). 

The  point  initiated  case  is  shown  in  Figure  2.  The  detonation 
wave  forms  a circular  arc  and  there  is  no  real  sign  of  a reflected 
wave.  This  may  be  explained  as  follows.  In  Figure  5,  the  detonation 
originates  at  the  point  P and  reaches  the  apex  of  the  conical  liner 
at  the  point,  0.  At  some  later  time,  the  detonation  front  will  be  a 
circular  arc  centered  at  P and  intersecting  the  wall  at  the  point,  Q. 
The  front  will  be  perpendicular  to  the  radius,  PQ,  and  make  an  angle, 
ili,,  with  the  wall.  The  base  angle  of  the  cone  is  a = 21°.  From  the 
law  of  sines: 


And , 


T esc  4> 
sin  <)> 

from  the  figure. 


= D*t  esc  (l80-a) 

= T sina/D#t 

cos  = sin  ( 90-^  ) = sin 


cos  iJj ^ = T sina/D*t 


Using  the  values:  D (detonation  velocity)  = 0.736  cm/psec, 

t (time)  = 10  psec,  T (see  Figure  5)  = 2.563  cm,  and  a = 21°  we 

see  that  must  be  82.8°.  Using  this  value  of  ^ and  the  detonation 

properties  of  pentolite,  the  coefficients  of  the  cubic  equation  of 
Appendix  A can  be  calculated,  and  the  resulting  equation  has  no 
positive  real  roots.  Consequently,  there  is  no  solution  to  the 
regular  reflection  problem.  The  reflected  shock  must  be  a Mach 
reflection  which  could  easily  be  obscured  by  the  relief  waves  from 
the  sides  (the  low  pressure  regions  on  the  right  of  Figure  2). 

In  the  plane  wave  initiated  case  (see  Figure  3)  the  detonation 
front  is  a straight,  horizontal  line  except  near  the  edges  of  the 
explosive  where  the  relief  wave  makes  itself  felt.  The  reflected 
wave  is  found  by  moving  horizontally  across  the  figure.  As  we 
move  from  right  to  left,  the  pressure  rises  and  then  falls  off. 

The  pressure  maxima  (found  in  the  shaded  boxes  of  Figure  3) 
approximately  locate  the  reflected  shock.  Near  the  liner,  the 
reflected  shock  makes  an  angle  of  nearly  45°  with  the  detonation 
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front  and  then  turns  to  a nearly  vertical  position  (90°  to  the 
front) . The  angle  of  incidence  is  90°-21°  = 69°  and  the  calculation 
of  Appendix  A shows  that  there  are  no  solutions  to  the  regular 
reflection  problem.  The  shock  system  we  are  looking  at  is  a Mach 
configuration  somewhat  obscured  by  the  smearing  out  of  the  detonation 
wave  over  three  zones. 

In  the  circumferentially  initiated  case  (Figure  4) , the  shock 
is  incident  upon  the  liner  head-on,,  and  the  reflection  corresponds 
to  a normal  reflection  fairly  well. 

Having  studied  the  shock  wave  systems,  we  must  now  turn  our 
attention  to  the  effects  of  these  shocks  upon  the  formation  of  the 
shaped  charge  jets.  The  first  and  easiest  comparison  between  the 
three  modes  of  initiation  is  the  speed  with  which  they  move  the 
metal.  This  can  be  summarized  by  looking  at  the  total  kinetic 
energy  transmitted  to  the  liner  (Table  1)  or  the  total  momentum 
(Table  2)  . 

Table  1 


Total  Kinetic  Energy  of  the  Metal  (ergs) 


10  ysec 

20  ysec 

30  ysec 

40  ysec 

Point 

7.642X1011 

7 . 132xl012 

2. 820xl013 

3 . 576xl013 

Plane 

2. 131xl012 

1. 623xl013 

4. 178x10 13 

4.246xl013 

Circum- 

ferential 

3.507X1011 

1.725xl013 

3.924xl013 

4. 142xl013 

Table 

2 

Total  Momentum  of  the 

Metal (gm  cm/sec) 

10  ysec 

20  ysec 

30  ysec 

40  ysec 

Point 

8.699xl06 

5. 098xl07 

1.469xl08 

1.907xl0S 

Plane 

1.617xl07 

9 . OlOxlO7 

1.980xl08 

2.306xl0S 

Circum- 

ferential 

3.009xl06 

8. 890xl07 

1. 888xlOS 

2.228xlOS 

These  values  are  plotted  in  Figures  6 and  7.  The  most  obvious 
result  is  the  inferiority  of  point  detonation  as  a means  of 
initiating  shaped  charges.  The  other  two  modes  of  initiation  are 
roughly  comparable,  with  the  plane  initiation  having  a slight 
advantage.  To  differentiate  further,  we  turn  our  attention  to 
the  jet  formed.  Figures  8,  9,  and  10  show  the  position  of  the 
metal  after  40  nsec.  As  can  be  seen  from  Figures  8,  9,  and  10, 
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in  all  three  cases,  at  40  yseconds,  the  liner  is  divided  into 
three  parts--a  slug  at  the  top,  the  undeformed  liner  on  the  right 
(which  is  approximately  in  its  original  shape)  and  the  jet  at  the 
bottom  which  is  separating  from  the  other  two.  In  each  figure  a 
horizontal  line  may  be  drawn  which  separates  the  jet  below  from 
the  slug  above.  This  line  is  drawn  where  the  surface  of  the  metal 
turns  sharply  to  the  side  and  begins  to  move  back  downward.  All 
the  cells  in  the  grid  below  this  line  belong  to  the  jet,  and  all 
those  above  do  not.  Taking  the  values  of  mass,  itk  ^ , and  velocity, 

v. .,  for  all  such  cells  and  forming  the  sums:  Em. Em. . v. .,  and 
13  2 1-3  1-3  13 

E 1/2  m^j  v^j  , gives  the  total  mass,  momentum,  and  energy  in  the 

jet.  The  results  for  the  three  different  modes  of  initiation  are 
summarized  in  Table  3 for  the  jet  at  30  yseconds,  and  in  Table  4 
for  the  jet  at  40  yseconds. 


Table  3 


Properties 

of  the  Jet  (30 

ysec. ) 

Jet  Mass 

Momentum 

Kinetic  Energy 

gm 

gm  cm/ sec 

ergs 

Point 

6.028 

0 , ?845xl07 

1. 228xl012 

Plane 

15.570 

1 . 111x10  7 

3.996xl012 

Circum- 

ferential 

16.458 

1 . 292xl07 

5.117xl012 

In  Tables  1 and  2 it  is  seen  that  circumferential  initiation  and 
plane  wave  initiation  are  approximately  equal  in  terms  of  the 
momentum  and  kinetic  energy  given  the  liner.  In  Table  3,  it  is 
seen  that  circumferential  initiation  puts  more  mass,  momentum,  and 
energy  in  the  jet.  Since  the  jet  is  the  penetrating  mechanism, 
circumferential  initiation  is  clearly  superior.  The  superiority 
is  even  more  evident  at  40  yseconds  as  seen  in  Table  4: 

Table  4 

Properties 

of  the  Jet  (40 

ysec. ) 

Jet  Mass 

Momentum 

Kinetic  Energy 

gm 

gm  cm/ sec 

ergs 

Point 

14.91 

1.103xl07 

4 . 085xl012 

Plane 

29.23 

2 . 144xl07 

7.881xl012 

Circum- 

ferential 

34.52 

2. 53^xl07 

9. 380xl012 
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The  values  at  40  yseconds  were  used  to  test  the  sensitivity  of 
the  results  to  the  exact  location  of  the  horizontal  line  which 
separates  the  jet  from  the  rest  of  the  metal.  Accordingly,  the 
sums  used  to  define  jet  mass,  momentum,  and  kinetic  energy  were 
extended  to  include  one  extra  row  of  zones.  The  revised  figures 
are  in  Table  5. 


Table  5 


Properties  of  the  Jet  at  40  yseconds — Extra  Zones  Included 


Jet  Mass 

Momentum 

Kinetic  Energy 

gm 

gm  cm/sec 

ergs 

Point 

*7 

5. 571xl012 

Detonation 

22.89 

1.587x10 

Plane  Wave 

41.43 

2 . 84  8x10  7 

9 . 9 35xl012 

Circum- 

11. 20xl012 

ferential 

45.69 

3.171x10 

If  the  second 

row  of  Table  4 

was  divided  by 

the  first,  the  result 

is  about  1.9 

for  each  column 

indicating  a 90%  improvement  over  the 

point  detonation  mode.  The  same  operation  on  Table  5 gives  1.8. 

The  difference  is  not  believed  to  be  significant.  For  the  quotient 
of  row  3 by  row  2,  Table  4 gives  1.2  and  Table  5 gives  1.1.  Again, 
the  difference  is  not  believed  to  be  significant  and  the  results  do 
not  seem  to  depend  critically  on  the  location  of  the  horizontal 
line  mentioned  above.  Even  more  significantly,  the  near  constancy 
of  the  ratio  seems  to  indicate  that  an  increase  in  mass  of  the  jet 
increases  the  other  quantities  proportionally.  The  distribution  of 
pressures  and  velocities  in  the  jet  at  30  yseconds  is  shown  in 
Figure  11  for  the  point  initiated  charge,  the  plane  wave  initiated 
charge,  and  the  circumferential  initiated  charge.  The  circum- 
ferentially initiated  charge  has  produced  a fully  developed  jet, 
because,  as  we  move  to  the  right  in  Figure  11  (toward  the  jet  tip) , 
the  velocity  increases  monotonically  and  the  pressure  falls  to  zero 
almost  immediately.  The  plane  wave  initiated  jet  is  not  yet 
developed  because  the  velocity,  though  increasing  rapidly  at  first, 
peaks  and  declines  somewhat.  The  pressures  (shown  below)  are  still 
substantial  and  further  accelerations  will  take  place.  The  point 
initiated  charge  produces  the  worst  jet  of  all  with  the  velocities 
peaking  and  declining  sharply  while  the  pressures  are  still  high. 
This  jet  obviously  needs  more  time  to  develop. 

The  situation  at  40  yseconds  is  shown  in  Figure  12.  The 
circumferentially  initiated  jet  has  not  changed  qualitatively. 

The  velocity  is  still  increasing  as  we  move  out  on  the  jet  and  the 
pressure  falls  to  zero  almost  immediately.  The  plane  wave  initiated 
jet  shows  signs  of  further  development.  The  velocity  no  longer 
peaks  and  falls  off  sharply,  but  falls  off  more  gradually.  Although 
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the  velocity  peak  has  risen,  the  tip  velocity  is  within  98.11% 
of  the  p ak  velocity  at  40  yseconds  as  opposed  to  only  96.64 °L 
at  30  yseconds.  The  pressures  also  have  fallen  off  sharply.  The 
point  initiated  jet  seems  to  have  undergone  very  little  further 
development.  The  velocity  gradient  is  still  wrong  and  the  pressures 
remain  high  even  up  to  the  tip  of  the  jet.  In  fact,  the  situation 
seems  to  be  getting  worse,  since  the  tip  velocity  is  further  away 
from  the  peak  velocity  at  40  yseconds  than  at  30  yseconds. 

From  these  results,  we  can  conclude  that,  for  the  charge 
configuration  considered  here,  the  circumferential  mode  of  initiation 
is  clearly  the  superior  mode.  The  point  detonation  mode,  although 
convenient,  is  a poor  way  to  detonate  this  charge.  (The  point 
initiated  case  is  not  discussed  in  the  following) . 

The  stability  of  the  remaining  two  cases  is  now  examined. 

It  is  known  that  a jet  will  be  radially  unstable  and  show  early 
signs  of  breakup  if  the  velocity  with  which  it  enters  the  stagnation 
point  is  supersonic2.  The  stagnation  point  is  easily  located 
since  it  is  in  the  zone  (on  the  axis)  with  the  highest  pressure 
(see  Figure  13) . The  sound  speed  at  this  point  can  ho  computed 
from  the  (Tillotson)  equation-of-state  by  the  relation 


c2  = 


3f 

9P 


3f 

31 


P/P: 


where  f is  the  Tillotson  function1 , I is  the  internal  energy,  p the 
density,  p the  pressure  and  c the  sound  speed.  The  difference  in 
velocity  between  the  jet  tip  and  the  velocity  at  this  point  (Vo)  must 
be  less  than  the  sound  speed.  For  the  plane  wave  initiated  charge, 
the  sound  speed  at  the  stagnation  point  is  found  to  be  .80  cm/ysec. 
The  jet  emerges  with  a velocity  of  .83  cm/ysec  and,  using  the 
Birkhoff,  Taylor,  Pugh,  and  MacDougall  theory3,  the  collapse 
velocity,  VQ  (see  figure  13),  can  be  found  from: 


V = V csc(B/2)  cos  (a/2) 

j ° 


Solving  the  second  of  these  equations  for  8,  and  eliminating 
6 from  the  first,  yields 


V 

o cos 
Id 


V 

^ cos2  (a/2)-sin  (a/2) 


(a/2) 


1/2 
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Solving  for  Vq : 


sin  (a/ 2) 


V = 
o 


cos  a cos2  (a/2) 
70 VT 


+ 


sin2  a 


4V 


•i/2 


Now  a = 21°,  Vj  - 0.83  cm/ysec,  and  d=0,736  cm/ysec,  so  that 

Vq  = 0.32  cm/ysec,  which  is  definitely  less  than  the  sound  speed 

computed  for  the  stagnation  zone.  For  the  circumferentially 
initiated  charge,  the  jet  tip  velocity,  V. , is  1.0  cm/ysecond. 

The  calculated  sound  velocity  at  the  stagnation  reqion  is  0.80  cm/ 
ysec  and,  from  the  above  formula,  V =0.48  cm/ysec.  This  case 
is  also  stable. 

We  can  conclude  that  circumferential  initiation  is  a 
superior  mode  of  initiation  producing  a better  and  faster  jet 
than  either  of  the  other  two  modes  of  initiation. 

PART  II  - VARIATION  OF  EXPLOSIVE  PARAMETERS 


The  question  of  the  effects  of  initiation  on  the  jet  produced 
is  part  of  a broader  question  of  the  effects  of  different  explosives 
on  jet  production.  More  specifically,  the  explosive  can  be 
characterized  by  its  loading  density,  pQ ; and  its  chemical  energy, 

Q (or  Ep . , the  energy  of  the  explosion  products  at  the  Chapman- 
Jouguet  state)  , and  by  the  parameter,  y,  that  occurs  in  the 
equation-of-state  for  the  reaction  products: 


p = (y-1)  E/p. 
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A systematic  study  of  the  importance  of  these  three  parameters 
was  undertaken  to  determine  the  effects  of  varying  each  of  them  upon 
the  jet  produced.  The  plane  wave  initiated  charge  (of  Part  I) 
with  pentolite  explosive  was  chosen  as  a standard  for  comparison. 

This  will  be  referred  to  as  the  baseline  design.  Three  more 
explosives  were  obtained  by  increasing  each  of  the  parameters 
mentioned  above  by  10%.  They  are  called  the  density  varied  explosive 
(Ap),  the  energy  varied  explosive  (AE) , and  the  gamma  varied 
explosive  (Ay) . The  jets  that  each  of  these  explosives  produces 
are  named  for  the  corresponding  explosive. 

The  results  of  the  calculation  can  be  summarized,  first  with 
regard  to  the  total  energy  and  momentum  (communicated)  to  the  metal, 
by  Tables  6 and  7;  and  with  regard  to  the  jet  properties  by  Table  8. 


Table  6 

Total  Kinetic  Energy  of  the  Metal  vs.  Time 


Kinetic  Energy  (ergs) 


Time 

Baseline 

AE 

Ay 

A p 

10 

yseconds 

0. 2131xl013 

0.  2058xl013 

0. 266xl013 

0.  239xl013 

20 

yseconds 

1. 623xl013 

1. 887xl013 

1. 966xl013 

1. 783xl013 

30 

yseconds 

4.178xl013 

4. 775xl013 

4.787xl013 

4.514xl013 

40 

yseconds 

4 . 246xl013 

4 . 693xl013 

4. 539xl013 

4. 599xl013 

Table  7 

Total  Momentum  of  the  Metal 

vs.  Time 

Momentum  (gni  cm/sec) 

Time 

Baseline 

AE 

Ay 

A P 

10 

yseconds 

0. 1617xl08 

0.1768xlOS 

0.1847xl08 

0.1734xl0S 

20 

yseconds 

0.  9010xl08 

0.9919xl0S 

1.014xl0S 

0.9574xlOS 

30 

yseconds 

1. 980xl08 

2.150xl0S 

2.129xlOS 

2.091xl0S 

40 

yseconds 

2.306xl0S 

2.471xlOS 

2.410xl0S 

2.435xlOS 
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These  results  are  displayed  in  Figures  14  and  15.  All  three  of 
the  varied  explosives  produce  a greater  kinetic  energy  and  momentum 
in  the  liner  than  does  the  original  pentolite,  but  there  is  little 
to  choose  from  between  the  three  of  them.  Accordingly,  we  turn  to 
Table  8 to  compare  the  jets  produced: 


Table  8 


Properties  of 

the  Jet  at  40 

yseconds 

Mass 

Momentum 

Kinetic  Energy 

gm 

gm  cm/ sec 

ergs 

iseline 

29.23 

2. 144xl07 

7. 881xl012 

AE 

31.26 

2.372xl07 

9 . 024xl012 

Ay 

30.43 

2.336xl07 

8 . 995xl012 

Ap 

36.63 

2.683xl07 

9. 863xl012 

From  Table  8 it  is  clear  that  an  increase  of  density  is  far  and 
away  the  most  effective  way  to  increase  jet  performance.  An  increase 
in  explosive  energy  is  next  and  an  increase  in  gamma  is  last, 
although  not  very  far  behind.  The  pressures  and  velocities  at 
30  yseconds  are  plotted  in  Figure  16  and  the  values  at  40  yseconds 
are  plotted  in  Figure  17.  In  these  figures  it  is  easily  seen  that  the 
four  curves  are  qualitatively  the  same.  But  there  seems  to  be  a 
discrepancy  between  Figure  16  and  Table  8 because  the  table  asserts 
that  the  density  variation  gives  the  highest  momentum  and  kinetic 
energy  while  the  figure  shows  that  the  density  variation  yields  a 
velocity  curve  that  lies  below  the  other  two.  In  fact,  from  the 
table,  the  order  of  the  jets  would  be  density  variation,  energy 
variation,  and  gamma  variation  while  the  figure  shows  just  the 
opposite:  gamma,  energy,  and  density.  The  key  to  this  puzzle  lies  in 

the  mass  going  into  the  jet.  The  reason  for  the  density  variation 
jet  being  so  much  higher  in  energy  and  momentum  than  the  gamma 
variation  jet  in  the  table  is  because  the  former  is  much  more  massive 
and  hence,  everything  that  depends  on  mass  is  much  higher  and  the 
velocity  difference  is  not  sufficient  to  offset  the  mass  difference. 

As  an  application  of  all  this,  consider  some  quantity,  J,  which 
characterizes  the  jet.  It  is  desired  to  estimate  the  change  in  J 
if  some  other  explosive  were  substituted  for  pentolite.  First 
calculate : 


C 

P 


J (pQ  + 0.1pQ)  - J (pQ) 
O.lpo 


and  then  calculate  Cy  and  by  similar  formulas, 


Then,  approximately 
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AJ  = Cp  Apo  + Cy  Ay  + CQ  AQ 

This  will  be  an  estimate  of  the  change  in  J when  pq,  y,  and  Q are 

changed  by  ApQ,  Ay,  and  AQ,  respectively.  As  an  example,  we  take 

the  experimentally  measured  jet  tip  velocities  measured  by  Simon^ 
for  the  BRL  precision  shaped  charge.  First,  we  need  the  calculated 
jet  tip  velocities  using  the  four  simulated  explosives: 


Table  9 

Jet  Tip  Velocity  vs.  Time  for  All  Four  Explosives 


Tip  Velocity  (cm/usec) 


Time 

Pentolite 

Ap 

AE 

Ay 

10 

. 36822 

.42040 

.41063 

.39724 

20 

. 76304 

.77880 

.79239 

. 79404 

30 

. 81856 

. 83708 

. 85218 

.85513 

40 

. 83242 

. 84723 

.86376 

.87004 

These 

values  are 

plotted  in  Figure 

1 8 and  the 

curves  seem  co 

level  off  near  40  yseconds.  Hence,  the  40  ysecond  value  is  chosen 
as  the  tip  velocity.  We  can  now  calculate  our  coefficients: 


r - v (y+.ly)  ~ v(y) 
Y .ly 


.13532 


Similarly,  = 6.46186  and  Cp  = .08975.  The  results  of  Simon  are 

summarized  in  Table  10.  His  values  of  detonation  pressure,  detonation 
velocity,  and  loading  density  (columns  2,  3,  and  4 of  Table  10)  lead 
to  different  values  of  gamma  and  energy  (columns  5 and  6 of  Table  10) 
for  pentolite  than  those  used  in  the  calculations.  But  his  values 
are  used  for  consistency's  sake.  Column  1 contains  the  name  of  the 
explosive  used;  columns  2,  3,  and  4 contain  estimated  values  of  the 
detonation  pressure,  detonation  velocity,  and  density  supplied  by  the 
installation  that  loaded  the  explosives.  Columns  5 and  6 contain 
the  values  of  y and  Q calculated  from  columns  2,  3,  and  4.  Columns  7, 
8,  and  9 contain  the  difference  between  the  values  of  pQ,  y»  and  Q 
for  the  explosive  of  column  1 and  the  values  for  pentolite.  Column  9 
contains  the  change  in  value  of  jet  tip  velocity  from  pentolite  as 
calculated  by  the  method  outlined  above,  and  column  11  contains  the 
change  in  jet  tip  velocity  as  measured  by  Simon.  Column  12  has 
the  percentage  difference  between  columns  10  and  11.  Although  the 
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percentage  difference  may  be  as  great  as  38%,  some  of  them  are 
rather  small  and  only  once  does  a difference  greater  than  28%  appear. 
This  is  fairly  good  when  allowance  is  made  not  merely  for  the  usual 
experimental  errors  involved  in  measurement  of  so  complex  a 
phenomenon  but  also  for  the  fact  that  the  coefficients  were  calculated 
for  a substantially  different  liner  than  the  one  used  in  the 
experiments.  Furthermore,  slight  errors  in  columns  2,  3,  and  4 lead 
to  relatively  large  errors  in  columns  8,9,  and  10.  All  things 
considered,  the  agreement  is  remarkably  good. 

CONCLUSIONS 


We  can  conclude  that  the  method  of  estimating  variations  in  shaped 
charge  performance  by  using  the  calculated  coefficients  seems  to  work 
and  merits  further  study. 

We  can  also  conclude,  all  other  things  being  equal  (which  rarely 
happens) , that  the  best  shaped  charge  explosive  is  the  densest. 

We  can  also  conclude  that  point  initiation  is  a poor  initiation 
system  and  that  the  circumferential  initiation  seems  to  be  the  best 
and  merits  further  study. 
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AXIAL  DISTANCE  (CM) 
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FIG.  2 PRESSURE  MAP  FOR  POINT  INITIATED  MODC,  10,.SEC 
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FIG.  3 PRESSURE  MAP  FOR  PLANE  WAVE  INITIATED  MODE.  IOmSEC 
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PRESSURE  RANGES  (dynes/sq.  cm) 
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FIG.  4 PRESSURE  MAP  FOR  THE  CIRCUMFERENTIALLY  INITIATED  MODE,  IOjxSEC 


FIG.  5 ANGLE  OF  INCIDENCE  OF  THE  DETONATION  WAVE  ON  THE  LINER 
FOR  THE  POINT  DETONATED  MODE 
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DISTANCE  TO  BOTTOM  OF  GRID  (CM) — •»  DISTANCE  TO  BOTTOM  OF  GRID  (CM) — ► 

FIG.  11  VELOCITY  AND  PRESSURE  IN  THE  JET  AT  30fi  SEC  FIG.  12  VELOCITY  AND  PRESSURE  IN  THE  JET  AT  40^  SEC 

FOR  ALL  THREE  MODES  OF  INITIATION  (ALL  3 MODES  OF  INITIATION) 


STAGNATION  REGION 


FIG.  13  COLLAPSE  OF  THE  LINER 


The  detonation  wave  bends  the  liner  and  drives  it  inward.  The  velocity  with  which 
the  material  flows  into  the  stagnation  region  is  VQ>  the  collapse  velocity.  The  velocity 
with  which  it  flows  out  towards  the  jet  is  Vj  and  the  velocity  with  which  it  flows  out 
toward  the  slug  is  V$. 


VELOCITY  (CM/^iSEC) 
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FIG.  18  JET  TIP  VELOCITY  VS  TIME  (ALL  FOUR  EXPLOSIVES) 
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Appendix  A 


Suppose  a detonation  wave  is  incident  upon  a rigid  wall. 

A shock  wave  appears  at  the  intersection  of  the  wall  and  the 
detonation  wave.  This  shock  is  called  the  reflected  shock  wave. 

We  wish  to  determine  the  angle  between  the  reflected  shock  and 
the  wall.  The  geometry  is  as  shown  in  Figure  19.  The  angles, 
tjK  and  ipr  are  the  angles  of  incidence  and  reflection,  respectively. 

The  flow  is  from  left  to  right  and,  in  the  coordinate  system 
chosen  (with  the  point,  0,  at  rest) , the  material  moves  into  the 
detonation  wave  with  a velocity,  D csc(i^).  It  emerges  from  the 

detonation  wave  with  a velocity,  q^,  which  makes  an  angle,  3^, 

with  the  normal  to  the  detonation  front.*  The  pressure,  density  and 
energy  of  the  flow,  as  it  emerges  from  the  detonation  front,  are 
the  Chapman- Jouguet  values.  The  flow  is  incident  upon  the  reflected 
shock  at  an  angle,  a^,  to  the  normal.  We  seek  the  angle  of 

reflection,  ipr,  and  this  can  be  obtained  from  the  angle,  a^,  through 

the  following  reasoning: 

From  Figure  19, 

0)  + 3^  + = it 

a)  + tt/2  + C + tt/2  = 2 7T 

C + ^ + «r  = * 

Hence : 

^ + 'Pj.  + = it 

Therefore 


tan  iJj  = -tan(iK  + a1  +3,) 
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Using  the  work  of  Dunne  ^5^}  , 

tan  cot 

(formula  8,  op.  cit.)  and  hence  2 

(Y  + 1)  + Y tan  <|>. 

tan  (\p.  + 8.)  = 

(2y  + 1)  tan 


When  iJk  and  y are  given,  this  expression  is  known  and  the  tangent 

of  the  angle  of  reflection  can  be  computed  by  the  earlier  formula 
as  soon  as  a.  is  known.  To  compute  a.  we  employ  formula  23  of  Dunne's 
work  [5J  : 2 1 

Y cot  tan  + (y  + D cot  ci^  -tan 

cot  ot^  [y  tan3  + (y  + 1)  + cot  tan 


2 

(Y  + 1)  cos 

Y2  + (Y+  l)2  cot2  4^1 

L J 

Y - 1 = 0 

Y + 1 


This  is  an  equation  for  cot  0^  with  and  y as  parameters.  By 

clearing  the  fractions  and  grouping  powers  of  cot  together,  we 
obtain: 


cot3  ^2 y 2 + (Y  - 1)  (Y2  + (Y  + l)2  cot2  tan  ip.^ 

cot3  ^2(y  + l)2  cot3  \p^  (Y  tan3  ip^  + Y + 1)| 

+ cot  |(y  + l)3  cot3  + y2  (y  + 3)  | tan  \Jk 

+ 2y2  (Y  tan2  ipi  + Y + U =0 

This  is  a cubic  equation  for  cot  . Since  0 < < 1/2  it  , the 

possible  values  of  are  positive  real  numbers.  Examining  the 

coefficients  of  the  equation,  we  see  that  the  coefficient  of 
2 

cot  is  negative  and  all  the  other  coefficients  are  positive. 

By  Descartes'  Law  of  Signs,  the  number  of  positive  real  roots  is 
either  2 or  0.  It  is  easily  shown  that  there  is  exactly  1 negative 
real  root  and  we  are  left  with  the  well  known  result:  that  either 
there  are  two  possible  reflections  (the  strong  and  the  weak 
reflections)  or  there  is  no  possible  regular  reflection.  In  the 
latter  case,  the  roots  are  imaginary. 
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To  determine  the  nature  of  the  reflection  in  our  cases,  we 
consider  first  the  point  detonated  explosive.  Here,  iJk  = 82.8° 

and  y = 2.78.  Using  the  standard  algorithm  for  solving  cubic 
equations : 

1 1) . =8, 

l 

a = 
b = 
c = 
d = 

Z = y -b/3a;  z3  + 3Hz 

H = 

G = 

u3  = 1/2 (-G  +V  G2  + 4H3  ) = .010394482 
u = .2182399723 
v = -H/u  = -2.287361695 
Z1  = u + v = -2.06912173 

2Z2  = - ti  + v)  + (u  - v)^3*  i 
2Zj  = - (u  + v)  - (u  - vh/T  i 

Z2  and  Z^  are  complex  and  Z^  is  negative.  The  corresponding  values 

of  cot  a are  also  either  complex  or  negative  and  there  is  no  genuine 
solution  to  the  regular  reflection  problem.  The  reflection  must  be 
a Mach  reflection.  The  plane  wave  initiated  problem  is  solved 

similarly.  Here  = 48°  and  y = 2.78  (as  above). 


2.831 


y = 2.78 


235.4582114 
-81.15242593 
361.9410554 
2774.5337962 

+ G = 0 reduced  equation 

.4991937361 

11.95713568 


coefficients 


ij'i  = 48°  y = 2.78 

a = 55.34500733 
b = -167.0185289 
c = 98.24218914 

d = 111.4283686 


coefficients 
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Z = y - b/3a  Z3  + 3HZ  + G = 0 (reduced 

cubic) 

G = 1.763190504 
H = -.4201864861 

u3  = 1/2  (-G  +t/g^  + 4H3)  = -.  0431303117 
u = -.3506933513 

v = -H/u  = -1.198159259 
Z1  = u + v = -1.548852611 
2Z  2 = 1.548852611  + .8474659077i 
2Z 3 = 1.548852611  - .8474659077i 

As  in  the  previous  case,  all  the  roots  of  the  reduced  cubic  are  either 
negative  or  imaginary.  The  corresponding  values  of  cot  are  also 

negative  or  imaginary  and  there  is  no  genuine  solution  to  the  reflec- 
tion problem.  The  reflection  must  be  an  irregular,  Mach  reflection. 


34 


NSWC/WOL  TR  77-123 


DISTRIBUTION  LIST 


Commander 

Ballistics  Research  Laboratory 
Alberdeen  Proving  Ground,  MD  21005 
(Dr.  F.  I.  Grace) 

(Dr.  A.  Dietrich) 

(Mr.  J.  Harrison) 

(Mr.  B.  Scott) 


Commander 
ARRADCOM 
Dover,  NJ  °7801 

(Mr.  G.  Randers-Pehrson) 
(Dr.  J.  Pearson) 


Science,  Systems  and  Software  Corp. 
P.  0.  Box  1620 
La  Jolla,  CA  92037 
(Mr.  Sedgwick) 


Commander 

Naval  Sea  Systems  Command 
Washington,  D.  C.  20362 
(PMS406 ) 

( PMS406E) 

(PMS406E3) 

(SEA0333) 


Commander 

Naval  Ocean  Systems  Command 
San  Diego,  CA  92152 

(F.  Marshall,  Code  0609) 
(ALWT  Library  Code  4473-B) 


Army  Mobility  Equipment  Research  and  Development  Center 
Fort  Belvoir,  VA  22060 
(J.  W.  Bond) 


I 


35 


NSWC/WOL  TR  77-123 


Distribution  List  (Cont.) 


Commander 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22314  12 


NDWX$WC(W)-5405/l  (R«.  5-78) 


TO  AID  IN  UPDATING  THE  DISTRIBUTION  LIST 
FOR  NAVAL  SURFACE  WEAPONS  CENTER,  WHITE 
OAK  TECHNICAL  REPORTS  PLEASE  COMPLETE  THE 
FORM  BELOW: 

TO  ALL  HOLDERS  OF  NSWC/WOL/TR  77-123 

by  J.  P.  Coughlin,  Code  R-13 

DO  NOT  RETURN  THIS  FORM  IF  ALL  INFORMATION  IS  CURRENT 


A.  FACILITY  NAME  AND  ADDRESS  (OLD)  (Show  Zip  Code) 


NEW  AOORESS  (Show  Zip  Code) 


B.  ATTENTION  LINE  ADDRESSES: 


C. 


□ REMOVE  THIS  FACILITY  FROM  THE  DISTRIBUTION  LlST  FOR  TECHNICAL  REPORTS  ON  THIS  SUBJECT. 


r 


DEPARTMENT  OF  THE  NAVY 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING.  MD.  20910 


OFFICIAL  BUSINESS 
PENALTY  FOR  PRIVATE  USE,  $300 


POSTAGE  AND  FEES  PAID 
DEPARTMENT  OF  THE  NAVY 
DOD  316 
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